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Abstract 

We study the zeros in the complex plane of the partition function for the Ising model coupled 
to 2d quantum gravity for complex magnetic field and real temperature, and for complex tem- 
perature and real magnetic field, respectively. We compute the zeros by using the exact solution 
coming from a two matrix model and by Monte Carlo simulations of Ising spins on dynamical 
triangulations. We present evidence that the zeros form simple one-dimensional curves in the 
complex plane, and that the critical behaviour of the system is governed by the scaling of the 
distribution of the singularities near the critical point. Despite the small size of the systems 
studied, we can obtain a reasonable estimate of the (known) critical exponents. 
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1 Introduction 



The understanding of the phase transition of a statistical mechanical model has long been con- 
nected to the study of its partition function zeros. This line of research was pioneered by Yang 
and Lee and has subsequently been pursued by many authors (for recent work see and 
references therein). In the thermodynamic limit, the zeros of the partition function of such a 
system in the complex coupling space accumulate infinitely close to the critical coupling on the 
real axis. This defines, in the infinite volume limit, disconnected regions in the complex plane, 
with different analytic structure for the thermodynamic functions which describe the phases of 
the system. 

The critical behaviour in the neighbourhood of a continuous phase transition can be extracted 
from the behaviour of the density of zeros near its singular points ||2|, ^, ^. In the best-understood 
models, these points occur at the ends of lines on which the zeros accumulate in the infinite volume 
limit. Such singular points, even when they occur at non-physical values of the coupling, can 
be considered as ordinary critical points with distinctive critical exponents Q (see also |^). For 
a physical phase transition, these singularities coincide with the (real) critical coupling of the 
system. 

Substantial progress, also from a practical point of view, has been made by applying renor- 
malization group arguments to the motion of zeros, thereby deriving finite size scaling relations 
by which the critical exponents can be computed Q. This finite-size scaling technique applied to 
complex zeros has been proven powerful in numerical computations of critical exponents. 

A special class of statistical systems are spin systems defined on random surfaces. In addition 
to being interesting systems per se, they have received special attention for their relation to 
conformal matter systems coupled to two dimensional quantum gravity. At criticality, the spin 
systems are described by conformal field theories (CFT) coupled to gravity. The critical exponents 
of such conformal theories have been found by Knizhnik, Polyakov and Zamolodchikov (KPZ) 
to be simply related to the critical exponents of the same CFT defined in flat space Q. In 
particular, the critical Ising model coupled to 2d gravity can be recast in the form of a c = 1/2 
conformal field theory coupled to 2d gravity. In addition, the equivalence of the Ising model 
on a random surface to a hermitean two matrix model, discovered and studied by Kazakov and 
Boulatov |1C, 11|, renders it exactly solvable. The critical exponents were first computed this 
way, outside of the CFT context. The system is found to have a third order phase transition from 
a high temperature disordered phase to a low temperature ordered phase. The coupling between 
matter and the geometry of the surface is found to change the critical exponents associated to 
geometry only at the phase transition. 

The location of the partition function zeros of the Ising model in the complex fugacity plane, 
the so called Lee- Yang zeros, in the presence of gravity has been studied by Staudacher []l2| . 
The somewhat surprising discovery that the Lee- Yang zeros lie on the unit circle in the complex 
fugacity plane for each lattice size N will be confirmed below for slightly bigger lattices than 
those considered in [12|. As discussed by Staudacher, the Lee- Yang theorem Q guarantees that 
the roots of the partition function of the Ising model on a fixed lattice Zfiat(GAr) for a large class 
of not too pathological lattices Gn will lie on the unit circle, but it is not a priori expected that 
the roots of the partition function Z^, obtained by summing the Zflat(GAr)'s over all dynamical 
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random lattices Gn, 

Zn = X! ^^bx{Gn), (1-1) 

{Gn} 

will also be located on the unit circle. The summation in (1.1) corresponds to an integration 
over the additional quantum degree of freedom introduced by a fluctuating metric. The Lee- 
Yang theorem for this case remains unproven. A priori the zeros can be located on different 
curves or even on sets with more complicated topology, or on two-dimensional regions. The latter 
cases have been observed for the complex temperature zeros of 2d Ising and Potts models on 
hierarchical lattices [E^, for anisotropic lattices p^, and interestingly, also for a regular lattice 



made of different polygons |29|. (See also ref. where a very interesting, albeit somewhat 
inconclusive attempt at studying the partition function zeros of Ising spin glasses is made, and 
ref. 1 31 1, where the zeros for aperiodic systems were studied.) 

In this paper we use the same method as in in order to compute the exact partition 
function for systems with up to 14 squares. We compute the Lee- Yang zeros of the partition 
function and confirm that they are located on the unit circle of the complex fugacity plane. We 
also compute the zeros of the partition function for complex temperature (Fisher zeros) and find 
that they move on curves in the complex temperature plane as the lattice size is varied. This 
indicates that in the infinite volume limit, the zeros accumulate to form dense sets on these 
curves. We test the scaling relations derived by Itzykson, Pearson and Zuber in the context 
of regular lattices |^] and extract from there some (combinations of) critical exponents. These 
scaling relations are based on the assumption that the spin-spin correlation length diverges at 
the critical point. Despite the small size of the systems studied, we find reasonable agreement 



with the critical exponents given in 11 1. This result is expected from the non-trivial fact 



that the matter correlation length diverges at the phase transition [l^ (see |15] for an example 
where this is not true) . This way we provide strong evidence that the phase transition is governed 
by the singularities of the distribution densities aX H = Q and (3c of the Lee- Yang and Fisher 
zeros respectively as is the case for the Ising model on a fixed lattice. The pattern of zeros we 
observe confirms the observation in |14] that no antiferromagnetic transition occurs for the Ising 
model on a square lattice when the spins are placed on the faces of the lattice. Using duality, we 
can study the pattern of the Fisher zeros of the dual model where the spins are placed on the 
vertices. We observe that the set of Fisher zeros maps onto itself under (3 —(3. This indicates 
the existence of an antiferromagnetic transition whose critical exponents are identical to the ones 
of the ferromagnetic transition^. This confirms the results reported in [|l^]. 

One of the methods for computing complex zeros is provided by Monte Carlo simulations and 
histogramming techniques. These techniques have been applied successfully to systems defined 



in flat space |lq , 17 1. They have also been tried on several systems coupled to gravity (random 



surfaces coupled to extrinsic curvature. Potts models, 4d quantum gravity [^), but it proved 
difficult to observe the partition function complex zeros. Here, we apply the method to Ising 
spins on dynamical triangulations [|l^ ], and we are able to observe the zeros lying closest to the 
real axis for a range of lattice sizes with up to 256 vertices. In general, it was easier to find the 
Lee- Yang zeros than the Fisher zeros. In the latter case, we were able to locate only the first 

*We thank Desmond Johnston for pointing this out to us. 



2 



zero for any given lattice size. After determining the position of the zeros as a function of lattice 
size, we test their scaling and the critical exponents computed are in reasonable agreement with 
their exactly known values. 

In section § we define our model and review the expected scaling behaviour of the zeros on 
the assumption of a diverging matter correlation length. In section |^ we describe in detail the 
methods we use in order to compute the zeros, and in section ^ we test their finite size scaling 
behaviour. Section ^ contains some concluding remarks. 



2 The model 

We will study one Ising model on a dynamical square or triangular lattice. In the former case 
the spins are located on the N faces (squares) of the lattice, in the latter on the Ny vertices. The 
partition function for a fixed lattice Gat is given by 

Zfl,t(G^,/3,^) = ^e^S(->'^-'^^+^^«'^' , (2.2) 

and the partition function Z{N, (3, H) for the model coupled to quantum gravity is obtained by 
summing Z^^x_{GN-,(i-,H) over all lattices Gn with spherical topology and N faces 

Z{N,[i,H)= Za^t{GN,P,H). (2.3) 

{Gn} 

In eq. ( |2.2[) denotes a sum over the neighbouring pairs of faces (in case of a square lattice) 

or vertices (in case of a triangular lattice) of Gjy. J2{a} denotes the sum over all the possible 
spin configurations, (3 is the inverse temperature and H the magnetic field. J2{Gn} ^^^^ discrete 
analogue of summing over all possible metrics up to diffeomorphisms on a spherical surface of 
fixed volume, since Gn defines a metric on the discretized surface if we define each link to have 
the same fixed length. The class of lattices we consider include degenerate ones which have double 
links and vertices of order one. 

The system ( p.3| ) undergoes a third order phase transition whose universal properties are 
independent of the microscopic details of the lattice fl^, 11 1. For the square lattice, the critical 



inverse temperature is given by /5c = ln2 ~ 0.69314718, and for the triangular lattice, /3c = 
ilniji^ w 0.16030370 0. The critical exponents are given by Q 

P = 1/2 (2.4) 

7 = 2 (2.5) 

6 = 5 (2.6) 

udH = 3, (2.7) 

which are defined in the usual way by the behaviour of the magnetization M ~ |/? — /3c|^, 
M ~ H^^^ and the magnetic susceptibility xm ~ 1/3 — Pc\~"' in the critical region, is the 
spin-spin correlation length exponent ^ ~ \P — Pc\~'^ , and dn denotes the Hausdorff dimension 
of space. 
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For a given lattice Gn, the partition function ( p.2| ) can be written (up to a multiplicative 
constant) polynomial 

N Ni 

Zf,,t{GN,P, = E E Cmny"'c^ , (2.8) 
m=0 n=0 

where c = e~^^, y = e~^^, is the number of faces in Gn (when we put the spins on the faces) 
and A'^; the number of links, y is called the fugacity. In the following we will also use the notation 
u = and K = u — Uc- 

Eq. ( |2.8D is a polynomial in two variables, and its analytic structure is entirely determined 
by its zeros in the complex plane. If we fix the temperature /?, Zfii^t{Gj\f , (3, H) is a polynomial in 
y and its zeros were shown by Yang and Lee j^] to lie on the unit circle in the complex fugacity 
plane. We expect the theorem to be true for a large class of (not too pathological) lattices Gn- 
Then one can write the free energy of the system as 

N 

i^flat(G/v, /?, y) = -In n (y - , (2.9) 

k=0 

where a factor f3 has been absorbed into F^^x.{Gn ■, (3-,y) and an additive constant has been sup- 
pressed {y = e~'^^). Here ykiP) are the zeros of Z^^-i{Gn , (i, H) in the y plane, called the Lee- Yang 
zeros. In the thermodynamic limit these zeros form dense sets on lines, which are Stokes lines 
for -Fflat(Gooj /5, y)- All the relevant information is then encoded in the density of zeros pyl{(3, 0) 
[H = iO). The free energy per spin in the limit of infinite volume A'^ — > oo is then given by 

i^flat(G'oo,/3,y) = - r depYL{(i:0)\n{y-e-^''). (2.10) 

J — TT 

Here —26 denotes the polar angle on the unit circle of the partition function zero in the complex 
y-plane. This notation is motivated by the Lee- Yang theorem. The singularities of the partition 
function are expected to occur at the ends of the lines on which the partition function zeros 
condense. In the case of (p.2|), in the high temperature phase the zeros form a gap > such 
that pyl{P,G) = for |^| < 6*0. The points yo = e^^*^'' are the Lee- Yang edge singularities 
and they can be regarded as conventional critical points Q. The (only) characteristic critical 
exponent a is 

PYL{P,e)r^{e-eor, (2.11) 

which implies that M {9 — 9qY . In ||2^ it was shown that a = —1/6 in two dimensions. 

As we approach the critical temperature from above, the gap closes (^o ~^ 0) and the zeros 
pinch the real axis at y = 1. This signals the onset of the phase transition, since the partition 
function has different analytic behaviour in the two disconnected parts of the y-plane. The 
information on the universal scaling of the partition function is given by the scaling of pyl{P,G) 
near 6 = 0. At /? = /3c, 6*0 = and if we define M ~ H^/^ we obtain 

PYL{l3c,0)r^6'/^ . (2.12) 



4 



Other critical exponents can be derived from scaling relations discussed in and later in this 
section. For a finite system the zeros are never on the positive real axis since Cmn is positive and 
c > for real f3 in eq. ( |2.8| ). When the system is critical, i.e. when the correlation length ,^ ~ L 
(L is the linear size of the system), the zeros approach the real axis infinitely close as L — > oo 
and one can apply finite size scaling in order to extract critical exponents. 

If we fix y in eq. ( p.sp , we obtain a polynomial of order Ni in c. Its zeros in the c-plane are 
the so-called Fisher zeros For many systems the Fisher zeros are located on curves C, but 
this is not necessarily true in general. For the 2d (zero-field) Ising model on a square lattice, 
these curves C are two intersecting circles (both in the c-plane and in the tanh(/?) plane). One of 
the curves intersects the real positive axis at the physical critical point. Similarly to the case of 
Lee- Yang zeros, the Fisher zeros condense on these two curves in the thermodynamic limit and 
the free energy per spin is given, in this limit, by their density pf{P,H): 

F^^,{Goo,P,H) = - ^ d/3V(/5',^)ln(c(/3) -e-2/5'), (2.13) 

where (3' is a complex number. Near the critical temperature and at if = 0, 

p^(/3,if = 0)~|/3-/3e|-"+\ (2.14) 

which implies that F^^t{G oo, H = 0) ~ 1/3-/3^1"°+^. Therefore a is the specific heat exponent. 

The possibility of extracting scaling exponents from the study of complex zeros of the partition 
function relies on their scaling behaviour under the renormalization group |^]. Applying finite 
size scaling is a convenient method often used in numerical simulations. Simple scaling arguments 
1^] give the position of the jth Lee- Yang zero as 

jj2^2p5/{udH) ^ /^.(i^ATl/C^'rfil)) . (2.15) 

where we have substituted N^/'^h [qj- linear size L of the system (N is the volume and dn 
is the Hausdorff dimension), fj is an analytic function and the Lee- Yang theorem implies that 
fj{0) < 0. If we invert the above relation we obtain the positions of the Fisher zeros. Hence we 
can deduce that the jth Lee- Yang zero will scale as 

Hj ~ at-WWh) ^ (2.16) 

and the j'th Fisher zero as 

Kj ~ AT-i/C^-rfff) . (2.17) 

The gap of the Lee- Yang edge singularity will scale as 

hI^-CK"^^^, (2.18) 

where C > and Hq = iOq. Conversely, in the scaling region the dependence of the jth Fisher 
zero on the value of the (real) magnetic field will be 



) ■ (2.19) 
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Eq. ( |2.19| ) implies that the trajectories of the motion of the Fisher zeros in the K plane with 
varying real nonzero H will form an angle 

with the real K axis. Using eq. ( 2.15| ) a stronger scaling relation can be derived: 

H]{N/j f^/^^+^^ = F (i^(iV/jf /('^'^«)) , (2.21) 

where F is a universal analytic scaling function such that -F(O) < 0, Hj is the jth Lee- Yang zero. 
Similarly 

^J = [i^) ^"(0)' (2.22) 

where Kj is the jth Fisher zero. 

In jl^, we provided evidence that there exists a divergent correlation length associated to 
matter fields coupled to gravity in the range of central charge < c < 1. Based on these results, 
we expected that the observed zeros, for large enough lattices, might show a behaviour compatible 
with the scaling hypothesis. It is the main purpose of this work to provide evidence that this is 
indeed the case. 



3 Computation of Complex Zeros 

Our computation of Lee- Yang zeros (real /3, complex H) and Fisher zeros (real H, complex /3) 
of the partition function Z{N, (3, H) was done using exact results from matrix models and Monte 
Carlo multihistogramming techniques. We describe these methods in this section. 



3.1 Exact determination of partition functions 

The exact solution of the Ising model on a square dynamical lattice came from solving the planar 



limit of the two matrix model defined in [10, |ll[. It was noticed that in this limit the free energy 
of the model 



Fin,g,(3,H) 



In 



r (Pir (1)2 



exp 



-Tr [<Pl + 4- 2ccl)l4 + (5e^/n)0f + (ffe-^/n)^^)] ) (3.23) 



where (j)i^2 are nxn hermitean matrices, equals the grand canonical partition function Z{g, (3, H) 
for the Ising model coupled to gravity: 



oo 

hm —Fin, g, /?, H) = Z{g, ^,H)=J2 c^Z{N, /3, H) 



n— >00 fl 



(3.24) 



N=l 
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where c = —^gcj (1 — c ) . The solution for Zi^g^ /3, H) is given by pO|, pi 



Z(g,/3,i7) = iln 



where 



^ 9 3 



1 ^ ^ 2z(cosh H-\) 



(3.26) 



.(1-Z)2 (l-z2)2 

Using eq. (|3.25|) and eq. (|3.26|) we can expand Z{g^ (3, H) in powers of g and read off the coefficients 
Z{N,P,H) (corresponding to lattice size with N squares) in eq. (3.24). 



We have done this using the symbolic manipulation program Mathematica for N < 14. We 
obtain 

TV 4[f ] 

Z{N, P, H) = c-^'y-^ ^ ^ , (3.27) 

m=0 n=0 

where [.] denotes integer part. The same calculation for < 6 was done in ||l^ and our results 
are in complete agreement. Then the roots of Z(N,P,H) for either fixed /3 or fixed H were 
computed numerically. Note that, apart from a trivial multiplicative factor, the partition function 
is a polynomial and has a finite number of zeros in the complex c and y planes. 



As already discussed above, and in |12] for N < 6, the Lee- Yang zeros can be seen, somewhat 
surprisingly, to lie on the unit circle in the complex y-plane. No Lee- Yang theorem has so far been 
proven for the partition function of the dynamical lattice. Fig. |l| shows the results for lattice sizes 
8 < < 14, at the (bulk) critical temperature. This figure clearly shows that the zero closest 
to the real positive axis approaches the point y = 1 as we increase the lattice size, indicating 
a vanishing gap at /3c as — > oo. We will study this approach quantitatively using finite size 
scaling in the next section. 

In Fig. ^ we show the motion of Lee- Yang zeros with varying /3 for the = 14 lattice as we 
approach /?c from the hot phase. The zeros close up towards the point y = 1 {H = 0) on the real 
axis, reflecting the expected vanishing of the gap in pylW^O) for /? — > /?c in the infinite volume 
limit. 

We note from Fig. |^ and Fig. Q that also the Fisher zeros form curves in the complex plane. In 
Fig. ^a) we display the approach of the Fisher zeros for if = to the (physical and unphysical) 
critical points on the real axis, c = ±1/4 for increasing lattice size. There are also Fisher zeros on 
the imaginary axis, which flow to c = ibzoo. In Fig. Q(b) we show the same curves mapped onto 
the c-plane, where the tilde refers to the dual spin model. The usual duality relation is given by 
c = tanh(/3). We note that our model is not self-dual, so the critical point Cc = 1/4 is different 
from the dual critical point Cc = 3/5. Fig. |3|(b) suggests the existence of an antiferromagnetic 
phase transition for the Ising model on a square lattice with spins placed on the vertices W^- It 



happens at ct^"^ = 5/3. This corresponds to the unphysical value of c = — 1/4, reflecting the fact 
that no antiferromagnetic transition occurs for the Ising model coupled to 2d gravity on a square 
lattice with spins placed on the faces. We notice that the zeros in the c-plane are mapped exactly 
onto each other under the transformation c — > 1/c. This implies that the critical exponents of 
the antiferromagnetic transition are identical to the ones of the ferromagnetic transition, as was 
noted in |14|. 
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Fig. ^a) shows the flow of Fisher zeros with varying magnetic field H for a fixed lattice size 
N = 14. The zeros on the arcs flow away from the imaginary axis as H increases, while the 
zeros on the imaginary axis again appear to move towards c = itzoo (note that these points are 
mapped onto c = (—1,0) in Fig. |^(c)). For H = the zeros should pinch the real axis like in 
Fig. ^, but for H ^ they seem to avoid the real axis. This is compatible with the absence of 
a phase transition in the presence of a symmetry breaking fleld. In Fig. ^(c) we show the same 
flow in the c-plane. 



3.2 Multihistogramming of Monte Carlo data 

For systems with more than a few spins, the study of complex zeros must rely on numerical 
methods. Monte Carlo methods are very powerful in studying a system using the full hamiltonian. 
They can be used to study singularities of observables in the complex plane as was done in 
the pioneering work ref. [|l6[. With the development of (single) histogramming techniques the 
partition function could be calculated for a continuous region in the coupling space and its 
zeros determined in the complex plane |17|. Using multihistogramming [pO], where Monte Carlo 



data taken at different couplings are combined optimally, a more accurate determination of the 
partition function is possible over a wide range of couplings. This method for computing complex 



zeros was first used in [21|. Here we provide a brief description of the method, in order to describe 
the procedure we followed. For details see p^]. 

The Monte Carlo updating of the triangulations was performed by the so-called flip algorithm 
and the spins were updated using a standard cluster algorithm. One updating sweep of the lattice 
consisted of approximately A'^; accepted flips where A'^; is the number of links of the triangulated 
surface. After a sweep of the lattice we updated the spin system. All this is by now standard and 
we refer to |18, 1£] for details. We use the high quality random number generator RANLUX |34| 



whose excellent statistical properties are due to its close relation to the Kolmogorov K-system 
originally proposed by Savvidy et al. p3[| . 

The lattice sizes that we simulated ranged from 32 to 256 vertices. In order to minimize finite 
size effects, we took the spins to be on the vertices of the triangulations and included degenerate 



triangulations with double links and vertices of order one |19]. First we made a rough map 
of the partition function. Since the partition function is invariant under /?—>/? + ikTi/2 and 
H ^ H + z/c7r/2 it is not necessary to calculate \Z]\i\'^ for Im(/3) > 7r/2 or lm.{H) > 7r/2. In our 
calculation we took H to be purely imaginary. Then we scanned the region where we expected to 
see the first zeros, and ultimately we took a denser grid of measurements around the points where 
the partition function touched zero within errors. For the largest lattice, 55 million sweeps were 
performed at each of 11 values f3k around the location of the zeros, in order to get a sufficiently 
accurate determination of the partition function in this region. For each (3^ we computed energy 
and magnetization histograms h^{E,M), where E is the total energy and M the (absolute value 
of the) magnetization, from which one can obtain an approximate estimate of the density of 
states: 

p{E, M)=Y, Wk{E)pk{E, M)=Y, Wk{E)h''{E, M)e^'=^ . (3.28) 
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The coefficients Wk{E) are the appropriate weights for an optimal determination of p{E,M) 



""'^^^ El exp m - Pi)E + Fi) ' ^^'^^^ 

where n^'^ is the number of measurements at Pk,i, and t^^i the respective autocorrelation times. 
Fi is the free energy at /?/, determined self-consistently from 

e-. . Zm = E , . (3.30) 

Then the partition function Z{f3) and any observable {0{E, M)) ^ for real /3 is given by: 

Z{f3) = Y.PiE,M)e-^^ (3.31) 

{0{E,M))^ = J—Y,0{E,M)p{E,M)e-^''. (3.32) 

^^^^ E,M 

Lee- Yang zeros are computed from the minima of 



(cos(Im(F) M)) - i(sin(Im(/?) M)) (3.33) 



for real /3 and imaginary H and .^(/?) = Z{P,0). The Fisher zeros were computed from the 
minima of 

= (cos(Im(/?) i?)) - i(sin(Im(/3) i?)) (3.34) 

for complex values of (3 and H = 0. The errors in the partition function and in the position of 
the complex zeros were computed by a standard binning procedure. 

Table 1 and Table 2 contain our results for the zeros. We first checked that multihistogram- 
ming was working quite well with our data by looking at multihistograms of the specific heat and 
susceptibility. It was relatively easy to observe the first few Lee- Yang zeros. For Fisher zeros 
the computation was harder and in a similar way we could only clearly observe the zero closest 
to the real f3 axis. Delicate cancellations in |Zp between contributions from the two terms in eq. 
( 3.34 ) make the zeros located further away from the real axis disappear in the statistical noise. 
This sets limitations to the size of the surface on which we were able to observe Fisher zeros. 
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/3± A/3 


Im(iJ) ± AIm(i?) 


64 


0.1718 ± 0.0005 


0.203 ± 0.003 


96 


0.1666 ± 0.0004 


0.1484 ± 0.0002 




0.2114 ± 0.0004 


0.2454 ± 0.0004 


128 


0.1643 ± 0.0014 


0.1182 ± 0.0005 




0.2037 ± 0.0012 


0.195 ± 0.005 




0.2386 ± 0.0018 


0.260 ± 0.003 


256 


0.1596 ± 0.0010 


0.069 ± 0.003 




0.1850 ± 0.0017 


0.112 ± 0.004 




0.2069 ± 0.0009 


0.148 ± 0.005 




0.2262 ± 0.0009 


0.185 ± 0.001 




0.2497 ± 0.0034 


0.175 ± 0.005 



Table 1. Lee-Yang zeros observed using muUihistogramming. 





Re(/3) ± ARe(/?) 


Im(/3) ± AIm(/3) 


32 


0.1428 ± 0.0016 


0.1773 ± 0.0005 


64 


0.1481 ± 0.0049 


0.1348 ± 0.0022 


96 


0.1533 ± 0.0021 


0.1121 ± 0.0034 


128 


0.1556 ± 0.0300 


0.1030 ± 0.0050 


256 


0.1527 ± 0.0035 


0.0788 ± 0.0047 



Table 2. Fisher zeros observed using muUihistogramming. 



4 Verification of scaling relations 



In this section we discuss the extent to which the scahng relations eqs. ( p. 15 - 2^2^ ) hold for the 
Ising model on a dynamical lattice. There is no a priori reason to expect these relations to hold. 
As discussed earlier, for the fixed lattice, their validity is due to the divergence of the spin-spin 
correlation length in the critical region. For the Ising model on a dynamical lattice, we expect 
the correlation length associated to geometry fluctuations to diverge pi. 



^, 24 1, but the same 

need not necessarily be true for the spin-spin correlation length even at a continuous transition 
|15||. Recently however, numerical evidence indicated that this is indeed the case, and that the 



system behaves as an ordinary statistical system near a third order phase transition [13|. 

Numerical simulations (2^, 24, 13| further indicate that the Hausdorff dimension of the system 
is very close to 4. Prom this value we can estimate the linear size L = N^/'^'^ of the systems that 
we are studying and find that it is indeed quite small. In spite of substantial finite size effects, our 
results will provide evidence that the scaling relations are indeed satisfied by the motion of the 
complex zeros with varying lattice size, couplings, or the order of the zero. We can, in general, 
observe reasonable agreement with the known critical exponents (we will refer to these as the 
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KPZ exponents). In the cases where a deviation from the expected value is observed, we can 
observe an asymptotic approach to this value with increasing lattice size. 



4.1 Zeros of the exact partition functions for small square lattices 



At K = and large enough N eq. ( 2.16| ) implies 

P6 



ln\Hj 



In(iV) + Cj 



The slope on a log-log plot of Hj vs. N is expected to take the value 



135 



-5/6 



-0.8333... 



(4.35) 



(4.36) 



for all j. Fig. ^ and Table 3 show our results for the first three Lee- Yang zeros. The extracted 
exponent combination is in reasonable agreement with the KPZ exponents, especially for the first 
zero. The errors reported in Table 2 are not true statistical errors (which are meaningless in this 
case). They are computed from the standard formula giving least squares linear fit errors and in 
this case they are simply a measure of the systematic deviation of the points from a straight line. 
We will follow this practice for our fits throughout this section. 



j 


Slope 


1 


-0.871 ± 0.002 


2 


-0.935 ± 0.002 


3 


-0.951 ± 0.002 



Table 3. Predictions for the combination of critical exponents —-^^ from scaling of 
the first three Lee-Yang zeros of exactly known partition functions . 



Eq. (g) imphes 



InlK,- 



ln|nj 



Uc\ 



^lnfA)+C. 

vdu viv; 



Similarly, since \Kj\ oc — /3c| + 0{\(3j — /3c| ), we would expect that 



ln|/3,- 

where the slope should be given by 



He 



vdu 



udn 



-In 

1 
3 



(4.37) 



(4.38) 



(4.39) 



and the constant C should be independent of j. Table 4 shows the results of fits to eq. (|4.3^ ), 
and similar fits with replaced by A'^, the number of vertices of the lattice. The corresponding 
fits using the scaling variable K did not yield straight lines in a log-log plot, as can be seen 
from Fig. |^. In this figure, we have plotted \Ki\ (lower curve) resp. \/3i — /3c\ (upper curve) 
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vs. on a logarithmic scale. The difference between the scaling behaviour of the two different 
scaling variables (which are expected to be identical for very large systems), shows that we are 
not deep in the scaling region and that finite size effects are important for the system sizes that 
we consider. The same must be said about the difference between the slopes and intercepts C of 
Table 4. From eq. ( ^.381) it is expected that the slopes and C should take the same value for all j, 
so that the data points would fall on one and the same curve. The difference between the results 
obtained by using and Ny should asymptotically decrease as the lattice size goes to infinity. 
(In the table, n denotes the number of points included in the fit.) 



Quantity fitted 


j 


Slope 


C 


n 


ln|^ - /3e| vs. ln(i/iV) 


1 


0.327 ± 0.001 


0.295 ± 0.002 


8 




2 


0.377 ± 0.002 


0.490 ± 0.005 


4 


Ml3-f3c\ vs. Hj/K) 


1 


0.383 ± 0.002 


0.492 ± 0.006 


6 




2 


0.438 ± 0.004 


0.667 ± 0.009 


4 



Table 4. Slopes and intercepts from the fits to eq. { 4-3d[ ). The theoretical value for the 
slope is = 1/3. The dependence of the slope and of C on j indicates finite size effects, 
as does the difference between N and N^. n is the number of degrees of freedom in the fits. 



The scaling relation eq. ( 2.19| ) predicts the angle that the trajectories of the zeros will form 
with respect to the Ke{K) axis as H is varied for large L. In our case the angle is expected to be 



2P6 



36° . 



(4.40) 



This prediction is valid for large system sizes L and small magnetic field. In Fig. ^(a) we show 
the trajectories of the first Fisher zero for = 4, 6, 8, 14 and magnetic fields in the interval 
—0.2 < H < 0.2. As expected from eq. (|2.19| ), H = corresponds to the "turning point" closest 
to the Ke{K) axis. Fig. ^(b) and Table 5 show the corresponding values of tan^/^, calculated from 
fits in the small H range —0.03 < H < 0.02. Although ip does not reach its infinite volume value, 
Fig. 0(b) shows the approach to it with increasing lattice size. The breakdown of scaling for very 
strong H can be observed in Fig. ^(c) where a trajectory for \H\ between zero and 2.6 are shown. 



N 


tanip lb Atan^ 


6 


6.489 ± 0.162 


8 


3.295 ± 0.036 


10 


2.308 ± 0.001 


12 


1.792 ± 0.006 


14 


1.450 ± 0.002 


oo 


0.727 
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Table 5. Measured {N = 6, 8, 14) and expected (eq. ^.4(\) ) slopes of 
trajectories of zeros in the complex K plane as the magnetic field H is varied, 
as a function of the size of the lattice. 



We also checked the validity of the scaling relation given by eq. ( 2.21| ). Taking K real, we 



plotted Hj{N/j)^/^ vs. (/3 — (5c){N / jY/'^ . The scahng function F is shown in Fig. |8[ All zeros 
of order j > 1 for lattice sizes (10 < N < 14) are included in the same graph. We observe that 
there exists a range of /? in the hot phase where the data points lie approximately on a univeral 
curve defining the function F. 

By varying the value of the exponent —6/ {6 + 1) = —5/6 of the scaling variable A = j/N on 
the y-axis away from its KPZ value, we observed a broadening of the curve and in this way we 
obtained a (somewhat subjective) determination of the exponent combinations 

^ - 0.85 ±0.05, (4.41) 



5 + 1 vdn 

(the equality of these follow from the general exponent equalities and was used in the derivation 
of the scaling relation ( 2.21| )). This is in excellent agreement with the KPZ value of 5/6. 



4.2 Zeros determined from multihistogramming 

In this section we describe the scaling of the zeros determined numerically from Monte Carlo 
data. The lattice sizes investigated here range from 32 to 256 vertices (60-508 triangles). We 
were able to determine the first five Lee- Yang zeros for the biggest lattice size, while for the 
smaller ones at most three zeros were visible. For the Fisher zeros, only the zero closest to the 
real axis was observed for any lattice size. 

The Lee- Yang zeros we observed with multihistogramming are listed in Table 1. The least- 



squares fits to eq. ( 4.35 ) give for the first and second Lee- Yang zero, 

^ = 0.773 ± 0.013 (j = 1) (4.42) 

= 0.788 ± 0.033 (j = 2) . (4.43) 

udH ' ^ ' 

All data points with j <2 shown in Table 1 were included in the fits. Using instead of N for 
the volume of the system we obtain: 

[35 



vdH 

(55 



0.787 ± 0.013 (i = 1) (4.44) 
0.800 ± 0.034 (i = 2) , (4.45) 



vdH 

which gives a measure of the presence of finite size effects. The results are in quite good agreement 
with the expected value of f55/vdH = 5/6. From the exactly known value vdu = 3 we obtain 

(36 = 2.36 ± 0.04 (j = 1) (4.46) 
155 = 2.40 ± 0.10 (j = 2) . (4.47) 
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The exact value of this combination of exponents is 2.5. 

The value of udn extracted from eq. ( [4.37 ) is somewhat lower than 3, but still remarkably 



close to it given the small size of the systems considered. The plot of the observed zeros is shown 
in Fig. |lO| (ln|i^| resp. ln|/3 — /3c| vs. In(A^)). We also calculated the slope of ln|i^| vs. ln(A''^) 
and of ln|/3 — /3c| vs. ln(A^„). The results are shown in Table 6. We note that by discarding the 
smallest lattices the value of vdn approaches 3, albeit with an increasing error. For example, 
including only the three largest lattices we obtain (using N for volume): 

^ 0.350 ±0.067. (4.48) 



Quantity fitted 


l/{udH) 


ln\K\ vs. ln{N) 


0.392 ± 0.016 


ln\K\ vs. In(iV^) 


0.407 ± 0.017 


ln\(3-Pc\ vs. ln{N) 


0.386 ± 0.014 


ln|/3 - Pc\ vs. ln(iV,) 


0.401 ± 0.014 



Table 6. Finite size scaling for the first Fisher zero. 



5 Conclusion 

We have computed the positions of the singularities of the partition function for the two- 
dimensional Ising model coupled to gravity in the complex plane. Both Lee- Yang and Fisher 
zeros were studied. 



We verified the result in [12| that the Lee- Yang zeros are located on the unit circle of the 
complex fugacity plane. This presents us with the challenge of proving a corresponding Lee- 
Yang theorem for the case when a fluctuating metric contributes an additional quantum degree 
of freedom. We also observed that the Fisher zeros form one dimensional curves in the complex 
temperature plane. 

Given the small size of the systems we studied, we obtained reasonable agreement with scaling 
laws derived from ordinary renormalization group arguments. Although one cannot use our 
results to accurately determine the scaling exponents, the extracted exponent combinations show 



a reasonable agreement with their exactly known values |11]. This, together with the fact that 
a diverging matter correlation length exists in the system [13|, gives us confidence to conclude 
that the critical behaviour of the system is given by the scaling of the distribution of the complex 
zeros of the partition function. 
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Figure 1: Lee- Yang zeros in the complex fugacity {y = c ) plane for the Ising model on square 
dynamical lattices of varying size N = 8(+), 9(x), 10(*), ll(n), 12(b), 13(o), 14(») at the bulk 
critical temperature. 
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Figure 2: Lee- Yang zeros in the complex fugacity plane for a square lattice with = 14 at the 
critical point = e'^/^-^ = 0.25(+)) and in the hot phase (c = 0.35(x), 0.45(*), 0.55(n), 0.65( 
■))• 
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Figure 3: (a) Fisher zeros in the complex c plane (c = e^^^) for the Ising model with spins on 
the faces of square dynamical lattices of varying size {N = 8, 10, 12, 14) at = 0. The zeros 
move on arcs and on the imaginary axis. The zeros flow towards ibzoo on the Im(c) axis, and 
towards the critical points c = (1/4,0) and c = (—1/4,0) on the real axis as the thermodynamic 
limit is approached. The physical critical point is c = 1/4. (b) The trajectories in (a) shown 
in the complex tanh(/?) plane. Recalling the duality relation c = tanh(/?), this corresponds to 
Fig. 3(a) for the model with spins on the vertices (with c related to the dual inverse temperature 



/5 by c = e-^P). The zeros approach the points (—1,0), (3/5,0) and (5/3,0) as N 



oo. The 

point c = 3/5 corresponds to the ferromagnetic and the point c = 5/3 to the antiferromagnetic 
transition of the model. The zeros are mapped onto each other under 1/c. 
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Figure 4: (a) Fisher zeros in the complex c plane for a square dynamical lattice of size AT = 14 
and at varying magnetic field {y = e"'^^ = 1.0, 1.6, 2.5). For zero magnetic field, the zeros with 
Re(c) > end in the vicinity of the physical critical point c = (1/4,0). (b) The first quadrant in 
(a) shown magnified for several values of the fugacity y = 1.0, 1.6, 2.5, 5.0. (c) The trajectories 
in (a) shown in the complex c-plane for zero and nonzero magnetic field. The fugacity takes the 
values y = 1.0, 1.3, 1.6, 1.9, 2.5 and N = 14. 
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Figure 5: Scaling of the first three Lee- Yang zeros for the Ising model on square dynamical lattices 
of size N {j labels the zero). The slopes are expected to be given by the exponent combination 
—(36/{i'dH) = —5/6 {dn is the Hausdorff dimension). The fitted slopes are given in Table 3. 
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Figure 6: Illustration of the scaling relation \K\ oc A^^i/('^'^^^) for the first Fisher zero at H = 0. 
Due to finite size effects, the curve for \K\ is not a straight line. On the other hand, \P — Pd 
scales well and yields a slope of —0.327 ± 0.01. The expected slope is —1/3. 
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Figure 7: (a) Trajectories of the first Fisher zero in the complex K plane for nonzero magnetic field 
—0.2 < lm.{H) < 0.2. The trajectories are for small square lattices of size N = 4, 6, 8, 10, 12, 14 
in order from left to right. The angle ip that the trajectories form with the Re(i^) axis is expected 
to approach ijj = 7r/(2/35) = 36° as ^ oo. (b) tani/j versus N for small lattices. The dashed 
line corresponds to the expected asymptotic value of tanip for large N. (c) Breakdown of scaling 
for strong magnetic field (up to \H\ 2.6). The angle ^ changes drastically for strong H. 
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Figure 8: Check of the scahng relation for the jth Lee-Yang zero H'^\-'^^/^^+'^) = F(i^A~^/('"^«)), 
where A = j/N and F is a universal function. All values of j > 1 for lattice sizes in the range 
10 < A'^ < 14 have been plotted on the same graph. The points are expected to fall on the same 
universal curve for large j. 
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Figure 9: The first (j = 1) and second {j = 2) Lee- Yang zero observed using multihistogramming 
on dynamical triangular lattices in the size range 64 < iVy < 256 (124 < N < 508). The slopes 
of the lines should be given by —f35/{udH) = —5/6. 
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Figure 10: The first Fisher zero observed using multihistogramming on dynamical triangular 
lattices in the size range 32 < A^^, < 256 (60 < N < 508). The slopes of the lines are expected to 
be -l/iydn) = -1/3. 
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